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Abstract 

We determine the ratio Fr /F n in QCD with Nf = 2 + 1 flavors of sea quarks, based on a 
series of lattice calculations with three different lattice spacings, large volumes and a simu- 
lated pion mass reaching down to about 190 MeV. We obtain Fk/F^ = 1.192(7) s tat(6) sy st- 
This result is then used to give an updated value of the CKM matrix element | V ns \ . The 
unitarity relation for the first row of this matrix is found to be well observed. 



°Centre de Physique Theorique is "UMR 6207 du CNRS et des universites d'Aix-Marseille I, d'Aix-Marseille 
II et du Sud Toulon- Var, affiliee a la FRUMAM" . 
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1 Introduction 



An accurate determination of the Cabibbo-Kobayashi-Maskawa (CKM) matrix element | V^ s | is 
important, because it allows to put bounds on possible extensions of the Standard Model of 
Particle Physics which, in turn, are relevant to guide direct searches, e.g. those planned at the 
Large Hadron Collider (LHC) at CERN [|TJ. 

Since the kaon is the lightest particle with strangeness, it is not surprising that much of 
the recent progress in this field derives from precision studies of leptonic and semi-leptonic 
decays of kaons (see e.g. 0). The theoretical challenge is to link the experimentally observed 
branching fractions to fundamental parameters in the underlying theory, e.g. the CKM matrix 
elements in the Standard Model. In this step lattice QCD calculations can help by providing 
decay constants and transition form factors. For an overview on these activities see e.g. the 
summary talks on kaon physics at the last three lattice conferences 0, f|, 

In this paper we shall follow a proposal by Marciano || to derive |V^ S | from \V u d\, using a 
lattice determination of the ratio Fk / F^ of leptonic decay constants. More specifically, in 
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the l.h.s., even after dividing it by the radiative correction factor (the last bracket on the r.h.s.), 
is known to 0.4% precision, if I = /i is considered 0. Also Mk, M v and are known with a 
relative precision of 3-10 -5 , 3-10 -6 and 10~ 7 , respectively @. And |V^d| has been determined 
from super-allowed nuclear /3-decays with an accuracy better than 0.03% J7|, ||. Therefore, the 
limiting factor for a precise determination of |Vu S | via (|l|) is Fk/F^ - this ratio is typically 
determined with a precision of a few percent in present days lattice QCD studies. 

Here we present a state-of-the art determination of Fx/F n in QCD in which all sources 
of systematic uncertainty are properly taken into account. We have performed a series of 
dynamical lattice calculations with a degenerate up and down quark mass m u d and a separate 
strange quark mass m s , a scheme commonly referred to as Nf = 2 + 1. The strange quark 
mass is roughly held fixed close to its physical value (m s ~mP hys ), while the light quark mass is 
heavier than in the real world, but varied all the way down to values which make a controlled 
extrapolation to the physical mass possible. The spatial size L is chosen sufficiently large, such 
that the extracted Fr / F n values can be corrected, by means of Chiral Perturbation Theory, for 
(small) finite- volume effects. Since our calculation includes three lattice spacings, a combined fit 
to all simulation results yields a controlled extrapolation to the continuum and to the physical 
mass point. Details of our action and algorithm have been specified in ||, where evidence for 
excellent scaling properties is also given. 

This paper is organized as follows. Sec. 2 contains a survey of where our simulation points 
are located in the (M%, 2M] C — M%) plane, a discussion of which chiral extrapolation formulae 
might be appropriate, an account of our strategies to quantify (and correct for) cutoff and 
finite- volume effects, and a description of the overall fitting procedure that is used to determine 



Fr/F-k- In Sec. 3 the final result is given and compared to other unquenched calculations [10 
13 of F K /F n . In Sec. 4 our result is converted, by means of (|Ip, into a value for |Ki S |. We find 
that, within errors, the first row unitarity relation |Vud| 2 + | Vu S | 2 + |Kib| 2 = 1 is well observed. 
This, in turn, puts tight constraints on possible extensions of the Standard Model [ITl H. 
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Figure 1: Overview of our simulation points in terms of M n and y2M^ — M%. The former 
gives a measure of the isospin averaged up and down quark mass while the later determines 
the strange quark mass. The symbols refer to the three lattice spacings a~0.124fm ((3 = 3.3), 
a ~ 0.083 fm (f3 = 3.57), and a ~ 0.065 fm (/3 = 3.7), respectively, and the physical point 
is marked with a cross. Error bars are statistical only. The numbers given correspond to 
one of the 18 two-point function, time fit-intervals that we use in our estimate of systematic 
uncertainties (see Sec. for details). At (3 = 3.3 and (3 = 3.7 only a single value of the strange 
quark mass, close to the physical value, is considered. At (3 = 3.57, simulations are performed 
at three values of the strange quark mass (see text for details). 



2 Simulation and analysis details 
2.1 Simulation parameters 

Our gauge and fermion actions, as well as details of the algorithm that is used to simulate QCD 
with Nf = 2 + 1 dynamical flavors, have been specified in [§ . Here, it is sufficient to say that 
our action combines good scaling properties with a one-to-one matching of lattice-to-continuum 
flavor; this avoids the complications of a low-energy theory with unphysical flavor symmetry 
breaking. 

We adjust the quark masses m u d, m s and set the lattice spacing a using aM n , clMk and 
aMs- We adopt a mass independent scale setting scheme. This means that we extrapolate, 
for each coupling f3 = 3.3, 3.57, 3.7, the values aM n , aMx, aM^ to the point where any two of 
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~/3 am u d am s TFxT traj. aM^ clMk Fk/F^ 

-0.0960 -0.057 16 3 x32 10000 0.4115(6) 0.4749(6) 1.0474(5) 

-0.1100 -0.057 16 3 x32 1450 0.322(1) 0.422(1) 1.079(2) 

-0.1200 -0.057 16 3 x64 4500 0.2448(9) 0.3826(6) 1.113(3) 

' -0.1233 -0.057 24 3 x64 2000 0.2105(8) 0.3668(6) 1.137(6) 

-0.1233 -0.057 32 3 x64 1300 0.211(1) 0.3663(8) 1.130(5) 

-0.1265 -0.057 24 3 x 64 2100 0.169(1) 0.3500(7) 1.153(6) 



-0.0318 0,-0.010 24 3 x 64 1 650,1650 0.2214(7),0.2178(5 ) 0.2883(7),0.2657(5) 1.085(1), 1.057(1) 

-0.0380 0,-0.010 24 3 x 64 1350,1550 0.1837(7),0.1778(7) 0.2720(6),0.2469(6) 1.107(3), 1.092(2) 

-0.0440 0,-0.007 32 3 x 64 1000,1000 0.1348(7),0.1320(7) 0.2531(6),0.2362(7) 1.154(4), 1.132(4) 

-0.0483 0,-0.007 48 3 x 64 500,1000 0.0865(8),0.0811(5) 0.2401(8),0.2210(5) 1.22(1), 1.178(7) 



-0.007 32 3 x 96 1100 0.2130(4) 0.2275(4) 1.0223(4) 

-0.013 0.0 32 3 x 96 1450 0.1830(4) 0.2123(3) 1.0476(9) 

3.7 -0.020 0.0 32 3 x 96 2050 0.1399(3) 0.1920(3) 1.089(2) 

-0.022 0.0 32 3 x 96 1350 0.1273(5) 0.1882(4) 1.102(2) 

-0.025 0.0 40 3 x96 1450 0.1021(4) 0.1788(4) 1.137(6) 



Table 1: Parameters and selected results from our simulations. The errors quoted here are 
purely statistical. These results correspond to one of the 18 two-point function, time fit- 
intervals that we use in our estimate of systematic uncertainties (see Sec.^^ for details). In 
this particular analysis, the scales at (3 = 3.3, 3.57, 3.7 are a' 1 = 1616(20) MeV, 2425(27) MeV, 
3142(37) MeV, respectively. 



the ratios in M n /Mx/M^ agree with their experimental values. Alternatively we have used Q 
instead of H to set the scale and used the difference between the two methods as an estimate of 
the systematics of scale fixing (see Sec. |2.5| ). Our simulations do not account for isospin breaking 
and electromagnetic interactions; thus experimental inputs need to be corrected for these effects. 
We use Mf^ = 135 MeV, M£ hys = 495 MeV, M| hys = 1318 MeV and M£ hys = 1672 MeV with 
an error of a few MeV |lTJ . 

Using LO chiral perturbation theory as a guide, we fix the bare strange quark mass at a 
given j3, such that the computed 2M\ — M% approaches its physical value when extrapolated 
in Ml to the physical point. As Tab.0 indicates and Fig. || illustrates, only one bare strange 
quark mass is used at (3 = 3.3 and at (3 = 3.7, and it is very near the value required to reach the 
physical point. At (3 = 3.57, three values of the strange quark mass are considered to provide 
some lever arm for an interpolation to the physical m s . The fact that our simulation points 
lie above the physical value of 2M\ — M% is due in part to cutoff effects in the relation of 
this quantity to the vector Ward identity strange quark mass. Indeed, for (3 = 3.3 and 3.7, 
1M\ — extrapolates near its physical value at M„ = M% hys , and the slope with which it 
extrapolates tends to zero as (3 increases. 

Regarding m U( i, we cover pion masses all the way down to ~ 190 MeV. On every ensemble 
the ratio Fk/F^ is measured with the valence up/down or strange quark mass set equal to the 
corresponding sea quark mass, so only the unitary theory is considered. Results are collected 
in Tab. [I]. It is worth noting that this same dataset was successfully used to determine the light 



hadron spectrum ITS 



We now give details of how we extrapolate to the physical mass point, to the continuum, 
and to infinite volume. We continue with an explanation of our global fitting strategy and the 
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pertinent assessment of systematic errors. 



2.2 Extrapolation to the physical mass point 

From the discussion above it is clear that we need to extrapolate our results for F% / F n in m u d to 
m^ ys , while the strange quark mass is already close to mf ys . It is important to note that even 
the extrapolation in m u d is small, amounting to only a couple percent in Fx/F n . Thus, any 
reasonable, well-motivated functional form which fits our results with a good confidence level 
should give a reliable estimate of Fk/F^ at the physical point. In order to assess the theoretical 
error that arises in this extrapolation, we choose to invoke three frameworks to parameterize 
the quark mass dependence: 

(i) SU(3) chiral perturbation theory (xPT) [pO] , 
(m) heavy kaon 577(2) chiral perturbation theory [IB|, 



(m) "Taylor" extrapolations involving polynomial ansaetze ||. 

We now summarize the theoretical background of these functional forms. 

Ad (i): 577(3) assumes that the u, d and s quark masses are small compared to the 



scale of chiral symmetry breaking. At NLO, the ratio F^/F^ can be written [20] 



F K 1 ~ ^%{fM>g(f ) + jMjHog(fe) + \Ml-\Ml\ log(^^)} + ^MlL 
F - 1 - 3^%{2M^ log(f ) + Milog(^)} + j*MlL h 

where, for P = n,K, rj, 



(2) 



m p i ,Ml. 



LO 



and where Mp, at this order, can be taken to be the lattice measured masses, with M 7 
(4M|- — M%)/3, using the LO 577(3) relation. Note that in obtaining (0), we have used the 
freedom to reshuffle subleading terms between numerator and denominator to cancel the sea 
quark contributions common to Fr and F n , which are proportional to the low energy constant 
(LEC) L 4 . That leaves only one NLO LEC in (0), i.e. L 5 (fi), whose /i dependence cancels 
the one in the logarithms. Throughout, we shall keep \x fixed to 770 MeV, the value which is 
customary in phenomenology. 

In addition to Eq. (|J), we consider two more SU(3) expressions, which are equivalent 
to (0) at NLO: the one obtained by fully expanding (Q) to NLO and the one obtained by 
expanding the inverse of @ to NLO. The use of these three different forms is one of the ways 
that we have to estimate the possible contributions of higher order terms (see Sec. [2.5| ). 



A number of collaborations have reported difficulties and large corrections when fitting 
their results for F n and Fx to NLO SU(3) chiral perturbation theory (or its partially quenched 
descendent) around M 7r >400MeV [§j]. This suggests that mf ys might be the upper end of the 
quark mass range where NLO chiral perturbation theory applies. However, such a statement 
may depend sensitively on the observable and in this respect F^jF^ is rather special. It is an 
£77(3)-flavor breaking ratio in which the sea quark contributions cancel at NLO. Moreover, the 
NLO expressions fit our results for this ratio very well with only two parameters, F and L 5 ; 
the size of NLO corrections are perfectly reasonable, of order 20%; and the values of F and 
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Ls that we obtain are acceptable for a fit performed at this order []. In addition, omission of 
all SU(3) fits would shift our final result for Fx/F n by less than one quarter of our final error, 
and slightly reduce the systematic error. 

Ad (ii): Given the previous discussion, it is clear that the heavy kaon SU(2) framework 
of [|1^] is an interesting alternative, since it does not assume that m s /(A7cF ) <C 1 , but rather 
treats the strange quark as a heavy degree of freedom, which is "integrated out". Of course, 
while this may be a good approximation at the physical value of m u d, where m u d/m s <C 1, it 
may break down for larger values, when m u d<m s , as is the case for our more massive points. 



What is needed for our analysis is Eq. (53) of []TB for the pion mass dependence of Fk, 
together with the standard SU (2) formula for the mass dependence of F n 0] . In these expres- 
sions, every low-energy constant bears an implicit m s dependence which can be written as a 
power series expansion in m s —m^ hys about the physical point. In practice we are only sensitive 
to the m s dependence of the LO term. Thus, we can express the ratio as 

F - pi «°° 1 l-p^log(M) • < 4 > 

where F, F denote the two-flavor chiral limit of the pion and kaon decay constant, respectively, 
while A, A are the energy scales associated with the respective LECs Q. In (|J), /xqcd is a typical 
hadronic energy scale. 

This expression for Fx/F n has five m s -independent parameters (F, F, a//XQCD, A, A), four 
if the ratio is expanded and only NLO terms are kept (F, F, q;///qcd, A 3//4 /A 2 ). We find that 
the former leads to unstable fits and thus do not use it in our analysis. One of the ways that 
we use to estimate possible contributions of higher order terms (see Sec. |2.5|) , is to consider 
two SU(2) xPT expressions, which are equivalent to (f|) at NLO: the one, already discussed, 
obtained by fully expanding @j to NLO and the other obtained by expanding the inverse of 
(||) to NLO. For reasons similar to those discussed under (i), we leave the determination of the 
LECs, which appear at NLO in Eq. (|j), for future investigation. 

In order to use the parameterization (Q), one has to give a definition of m s -mf ys . SU(3) 
xPT, together with Fig. ^ and the discussion surrounding it, suggests that ([2M|- — M%] — 
[...]p hys )//iQCD oc (m s -mf ys )[l + 0(m s -mP hys )], up to negligible, physical, m ud and m^ ys 
corrections, in the range of quark masses covered in our simulations. Thus, we use this difference 
of meson masses squared as a measure of m s -mf ys . 

Ad (in): Both two and three flavor chiral perturbation theory are expansions about a 
singular point, (m ud ,m s ) = (0,0) in the case of SU(3) and (m u d,m s ) = (0,mP hys ) in the case 
of SU(2). However, here we are interested in the physical mass point, (mJ s ,mP hys ), which is 
nearer the region in which we have lattice results than it is to either singular point. Thus, it 
makes sense to consider an expansion about a regular point which encompasses both the lattice 
results and the physical point ||. As discussed above, the data at two lattice spacings (ft = 3.3 
and ft = 3.7) have been generated at a fixed value of am s , and the data at ft — 3.57 can easily 
be interpolated to a fixed am s . On the other hand, an extrapolation is needed in m u ^. It is 



1 Note that the ratio Fk/F^ is not well suited for a determination of Fq and L5. In this ratio, the LO LEC 
-Fb only appears at NLO where its value becomes highly correlated with that of L5. A serious determination of 
these constants would actually require a dedicated study which is beyond the scope of this paper. 

2 In the standard 51/ (2) framework the former is usually denoted A4, since it is associated with I4 |pifl. 
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natural to consider the following expansion parameters 



A. 
A* 



Ml - -(Mf 



cut\2 



/a*qcd 



M: 



K 



(Mt y r /a* 



''QCD 



(5) 
(6) 



with /iQCD as above and M^ ut the heaviest pion mass included in the fit. The discussion in 
|| suggests that the mass dependence of F K /F n in our ensembles and at the physical mass 
point can be described by a low order polynomial in these variables, leading to the "Taylor" or 
"flavor" ansatz 

^ = A + A,A n + A 2 Al + BiA K . (7) 

One of the ways that we use to estimate possible contributions of higher order terms (see 
Sec. |2.5| ), is to consider the same functional form for F^/Fk- Thus, we consider two flavor 
ansaetze. 

There are, of course, many possible variants of these ansaetze. For instance, Ak in Eq. (|6|) 
could be defined in terms of Mk instead of M K . One could also enforce SU(3) flavor symmetry, 
i.e. Fk/ 'F n = 1 when m u d = m s . This could be done, for instance, with an expansion of the 
form Fk/Fk = 1 + (M K — Ml) x [polynomial in A n and Ak]- However, the expansion of Eq. (0) 
provides a better description of our data than the alternatives that we have tried. 



2.3 Continuum limit 

To maximize the use of our data, we combine the continuum extrapolation with our interpola- 
tions and extrapolations to the physical strange and up-down quark mass point. As mentioned 
above, Fk/F^ is an S , f/(3)-flavor breaking ratio, so that cutoff effects must be proportional 
to m s —m ud . Although || suggests that they are quadratic in a, our action is formally only 
improved up to 0(aa s ) and we cannot exclude, a priori, the presence of linear discretization 
errors. Moreover, the effects are small enough that, despite a factor of almost two in lattice 
spacing, we cannot distinguish a from a 2 in our data. Thus, in our analysis, we consider 
three possibilities: no cutoff effects; cutoff effects of the form a(m s —m U( i), and cutoff effects 
proportional to a 2 iiQcv(m s —m U( i) . Again, /zqcd is a typical hadronic energy scale. 
In the case of SU(3) x?T, the cutoff effects can be accounted for by adding 



Fk 
F w 



a(M 2 K -Ml)/fi QC B 
cx I or (8) 
a\M 2 K -Ml) 



to the functional forms for the mass dependence of Fk/F^ given in the previous section. Here 
c is the relevant parameter. 

In the usual power counting scheme of SU(3) Wilson x?T, in which 0(a) = 0{m q ) p2[ , the 



cutoff effects considered in (^) are NNLO in the 0(a) case and even N 3 LO for those proportional 
to a 2 . This should be compared to the physical contributions to the deviation of Fk/F v from 
1, which are NLO. Thus, these cutoff effects are small parametrically and, as we will see below, 
they are also small numerically. In fact, they are consistent with zero within our statistical 
errors. Clearly then, any reasonable parameterization of these small S"?7(3)-flavor breaking, 
cutoff effects is sufficient for our purposes. Thus, we use the parametrizations of Eq. (Q) to 
subtract them, also in the context of SU(2) and flavor expansions. 



7 



2.4 Infinite volume limit 



The finite spatial size L of the box affects masses and decay constants of stable states in a 
manner proportional to exp(— M n L) f23| . In our simulations the bound M 7r L>4 is maintained. 
It turns out that the sign of the shift is the same for both decay constants so that such effects 
partly cancel in the ratio. In chiral perturbation theory to 1-loop order, F 7T (L)/F 7T has been 
calculated in [24] and Fk(L)/Fk in pB" |. At the 2- loop level, both ratios have been determined 
2B|. The generic formulae take the form 



in 



F K (L) 
F k 

Fir 



i + E 

n=l 
oo 

i + E 



m(n) 



Ml 



n 
m(n, 



M W L F K (4tt^) 2 
1 Ml 



n=l 



M W L (AnF n ) 



j(2) , 
1 F K + 



4? + 



M K rW | 



(4^ 



/i 4) + 



(9) 
(10) 



with m(n) tabulated in g§. With = — 4K 1 (y / nM 7r L) and = -^(^J^ii), where 



(2) 



ifi(-) is a Bessel function of the second kind, one obtains the compact 1-loop expression [P6j 



1 + 



3^ 



K 



E 

n=l 



min 



Ki(^/nM n L) 



(11) 



which is readily evaluated. On the other hand, the expressions It?,, •, I 



p^, 1$ in the 2- loop part are 
more cumbersome to deal with. Since finite-volume effects are independent of cutoff effects (see 
e.g. the discussion in |26[), we choose to first correct the data for finite volume effects before 
we actually perform the global fit (see below). The masses of the mesons are also corrected for 
finite volume effects according to p6| . 

To estimate the error associated with the finite volume effects, the difference between the 
1-loop shift of F K /F n and the corrected value of F K /F n using only the expression for F n 
(that should be an upper bound on the finite volume correction), is used as a measure of the 
uncertainty of the correction (see Sec. |2.5| ) . In our final analysis, the finite volume shift is 
smaller than the statistical error in most of our ensembles. 



2.5 Fitting strategy and treatment of theoretical errors 

Our goal is to obtain Fk/F^ at the physical point, in the continuum and in infinite volume. 
To this end we perform a global fit which simultaneously extrapolates or interpolates M% — > 
M2|p hys , M\ — > M|-| p h ys and a — > 0, after the data have been corrected for very small finite 
volume effects, using the two-loop chiral perturbation theory results discussed above. To assess 
the various systematic uncertainties associated with our analysis, we consider a large number 
of alternative procedures. 

To estimate the systematic uncertainty associated with setting the scale, we have repeated 
the analysis setting the scale with the S and the Q. 

To estimate the possible contributions of excited states to the correlators used, we repeat 
our analysis using a total of 18 different time intervals: t m i n /a = 5 or 6, for (3 = 3.3; 7, 8 or 9, 
for j3 = 3.57; 10, 11 or 12 for = 3.7. All of these intervals are chosen so that the correlation 
functions are strongly dominated by the ground state. However, the intervals which begin at 
earlier times may receive some small contamination from excited states. 
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Uncertainties associated with the necessary extrapolation to the physical up and down 
quark-mass point are assessed by varying the functional form used as well as by varying the 
range of pion masses considered. As discussed in Sec. pT2j we consider 3 forms based on the 
NLO SU(3) formula of Eq. (g), 2 forms based on the NLO SU(2) expression of Eq. (§, and 2 
flavor ansaetze, based on Eq. ([?]). Moreover, we impose 2 cuts on the pion mass: M m <350 MeV 
and 460 MeV. 

Our results for Fx/F n display a small dependence on lattice spacing. As described in 



Sec. |2.3| , to estimate the systematic associated with the continuum extrapolation we consider 
fits with and without 0(a 2 ) and 0(a) Symanzik factors. 

All of this amounts to performing 2-18-7-2-3 = 1512 alternative analyses. The central value 
obtained from each procedure is weighted with the quality of the (correlated) fit to construct 
a distribution. The median and the 16-th/84-th percentiles yield the final central value and 
the systematic error associated with possible excited state contributions, scale setting, and the 
chiral and continuum extrapolations. To obtain the final systematic error of our computation, 
we add a finite-volume uncertainty in quadrature to the error already obtained. This finite- 
volume uncertainty is given by the weighted (with the quality of the fit) standard deviation of 
our final central value and the ones obtained repeating the whole procedure with finite-volume 
effects subtracted at one loop and the upper bound computed as the 2 loop correction in F n 
only. We do not include these two alternative options in the set of procedures which yield our 
final central value because we know, a priori, that the two-loop expressions of |2B| are more 
accurate than the one-loop ones. To determine the statistical error, the whole procedure is 
bootstrapped (with 2000 samples) and the variance of the resulting medians is computed. 

There is a final source of theoretical error which we wish to comment on: the one associated 
with the fact that our calculation is performed with m u =m ( i and in the absence of electromag- 
netism. As discussed in the Introduction, we correct for these effects at leading order, up to 
a few per mil uncertainties. The latter have a negligible effect on F K /F W . One also expects 
direct isospin violation in Fk/F^. These were found to be negligible in fTI] . 



Having estimated the total systematic error, it is interesting to decompose it into its individ- 
ual contributions. By construction, the contribution from the uncertainty in the finite-volume 
corrections is already known. To quantify the other contributions, we construct a distribution 
for each of the possible alternative procedures corresponding to the source of theoretical error 
under investigation. These distributions are obtained by varying over all of the other procedures 
and weighing the results by the total fit quality. Then, we compute the medians and average 
fit qualities of these distributions. Clearly the spread of the medians measure the uncertainty 
associated with the specific source of error. We use the standard deviation of these medians 
weighted by the average fit quality as an estimate of the error under consideration. 

A "snapshot" fit (with a specific choice for the time intervals used in fitting the correlators, 
scale setting, pion mass range) can be seen in Fig. g To avoid the complications of a multi- 
dimensional plot, the extrapolation is shown as a function of the pion mass only. The data 
have been corrected for the deviation of the simulated m s from mf ys . In other words, what is 
shown is data(M2, 2M 2 K -Ml) -fit(M*, 2M^-M 2 ) + fit(M*, pMjH^W). Th e figure shows 
one flavor fit with no cutoff effects. We emphasize that y^/dof for our correlated fits are close 
to one. 

Fig.|3| shows the distribution of our 1512 alternative fitting procedures. As can been seen, 
the use of different formulas to extrapolate to the physical point only increases the systematic 
error. This is also true for the other sources of systematic error: scale setting, time intervals for 
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Figure 2: Extrapolation of the lattice data to the physical point for a particular choice of 
two-point function fits (t miQ /a = 6,8,11 for (3 = 3.3,3.57,3.7, respectively), mass cut (M n < 
460 MeV ) and using the S to set the scale. The plot shows one (of the 21 ) fits used to estimate 
the uncertainty associated with the functional form used for the mass extrapolation. The data 
have been slightly adjusted to the physical strange quark mass, as well as corrected for tiny 
Gnite-volume effects (see text for details). 



the fit to the correlators, pion mass ranges and cutoff effects. The figure also shows our final 
result with the total final error. 



3 Results and discussion 



Following the procedure outlined above, our final result is 



Fk 



phys 



1.192(7)^(6) 



syst 



or 



F K 



phys 



0.839(5) stat (4) syst 



(12) 



at the physical point, where all sources of systematic error have been included. Tab.0 gives 
a breakdown of the systematic error according to the various sources. We conclude that our 
main source of systematic error comes from the chiral extrapolation (functional form and pion 
mass range), followed by cut-off effects. 
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Figure 3: Distribution of final values for Fk/ 'F n . The large background distribution represents 
the values of Fk / F n obtained with different extrapolation formulas, pion mass cuts, parameter- 
ization of cutoff effects, time intervals and different methods to set the scale. Also shown is the 
final result (solid vertical line) and final systematic error interval (dashed vertical lines), which 
includes the finite volume error estimate. The smaller distributions group the values according 
to the chiral formula used for the extrapolation. The error associated with the chiral extrapo- 
lation is computed as the weighted (by total fit quality) standard deviations of the medians of 
the grouped distributions. 



In the same manner we may read off from our fits that the value in the 2-flavor chiral limit 



Fk 

Ft- 



m ud =0,m s phys 



1.217(10)^(9) 



syst 



or 



F K 



m ud =0,m s phys 



0.821(7)^(6) 



syst 



(13) 



differs by about 2% from the physical value, though with larger uncertainties. Upon combining 
this ratio in the SU(2) chiral limit with the phenomenological value F = lim mud _j.o F n = 
86.2(5) MeV 0, we obtain \im mud ^ F K = 104.9(1.3) MeV at fixed physical m s . 

In Fig.|] a survey of recent determinations of Fk/F v in unquenched QCD simulations is pre- 
sented. These include results by JLQCD [|Kj, MILC P, Q, NPLQCD |3|, HPQCD/UKQCD 
g, ETM RBC/UKQCD |6|, PACS-CS |7|, and Aubin et al. flTJJ. It is worth noting 



that these results show a good overall consistency when one excludes the outlier point of |L0 



The result by HPQCD/UKQCD remains the one with the highest claimed accuracy. Our 
precision is similar to that of another state-of-the-art calculation, the one by MILC. Given 
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Source of systematic error error on Fx / F n 
Chiral Extrapolation: 

- Functional form 3.3 x 10 -3 

- Pion mass range 3.0 x 10 -3 
Continuum extrapolation 3.3 x 10 -3 
Excited states 1.9 x 1CT 3 
Scale setting 1.0 x 10~ 3 
Finite volume 6.2 x 10 -4 

Table 2: Breakdown of the total systematic error on Fk/F^ into its various components, in 
order of decreasing importance. 



that we reach smaller pion masses, that more generally we explore a large range of simulation 
parameters which allows us to control all sources of systematic error, that we use an action 
with a one-to-one matching of lattice-to-continuum flavor and that we avoid the specifics of a 
partially quenched framework, our result solidifies the claim that Fx/F n is known to better 
than 1%. Moreover, agreement between such different approaches can only bolster confidence 
in the reliability of lattice calculations per se. 



4 Update on \V US \ and check of CKM unitarity relation 

With the result (|12"D in hand, we can now focus on CKM matrix elements. In this respect there 
are two options - we may assume the SM (and hence CKM unitarity) and determine |Vud| 
and |Vu S |, or we may use phenomenological input on |Vud| to derive |Vu S | and hence test CKM 
unitarity (under the assumption of quark-flavor universality) in a model-independent way. 

The first step, needed in either case, is to simplify Marciano's Eq. (|I]). The most recent 
update of the Flavianet kaon working group is [|J 

>!n>l / /v 0.27599(59) . (14) 



| Kid | Fit 

Combining this with our result ( |T2"D yields the ratio 

|Ks|/|Kd| = 0.2315(19) . (15) 
Now for the two options. If we assume unitarity, (15) and | Vub| = (3.93±0.36) 10~ 3 @ imply 



|K d | = 0.97422(40) , \VJ\ = 0.2256(17) . (16) 

On the other hand, if we combine (|5|) with the most precise information on the first CKM 
matrix element available today, jV^I — 0.97425(22) ||, we obtain (again) 

|K S | = 0.2256(18) . (17) 

Similarly, by also including the above mentioned result for \V U ^\ we find 

IKdl 2 + IKsl 2 + IKbl 2 = 1.0001(9) . (18) 
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Figure 4: Our result ( fE^ compared to previous unquenched lattice QCD calculations. The two 
error bars refer to the statistical and to the combined statistical and systematic uncertainties, 
respectively. The top two results were obtained with Nf = 2 simulations. The next four were 
computed on N f = 2 + 1 MILC conhgurations. RBC/UKQCD and PACS-CS determinations 
were obtained on distinct Nf = 2 + 1 ensembles. 



With the first-row unitarity relation (which is genuine to the CKM paradigm) being so well 
observed, there is no support for "beyond the Standard Model" physics contributions to these 
processes. Of course, with substantially improved precision on both the theoretical and the 
experimental side, this might change in the future. 



5 Summary 

Based on a series of 18 large scale lattice computations the ratio Fk/F^ has been determined 
in Nf = 2 + 1 flavor QCD at the physical mass point, in the continuum and in infinite volume. 
The overall precision attained is at the 1% level, with statistical and systematic errors being of 
similar magnitude. The systematic error, in turn, splits into uncertainties arising from the ex- 
trapolation to the physical pion mass, to the continuum, to possible excited state contributions 
and to finite volume effects. We find that the main source of systematic error comes from the 
extrapolation to the physical point. In this respect it is useful that our simulations cover the 
pion mass range down to about 190 MeV, with small cutoff effects and volumes large enough 
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to maintain the bound M 7r L>4. 

Following Marciano's suggestion ||, we use our result (|T^) for Fx/F n to obtain the value (|I5|) 
f° r |Kis|/|Kid| an d subsequently |U US |. In turn, these results are used to test first row unitarity, 
which we find is satisfied at the 1 per mil level, thereby imposing interesting contraints on "new 
physics" scenarios (see e.g. 0). 
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